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High-temperature superconductors exhibit a characteristic hourglass-shaped spectrum of magnetic 
fluctuations which most likely contribute to the pairing glue in the cuprates. Recent neutron scattering 
experiments in strongly underdoped compounds have revealed a significant low energy anisotropy of these 
fluctuations which we explain by a model in which topological defects of the antiferromagnet clump to 
producing domain wall segments with ferronematic order. This state does not invoke global charge order 
but breaks C4 rotational and inversion symmetry. The incommensurability of the low doping 
charge-disordered state is in good agreement with experiment and interpolates smoothly with the 
incommensurability of the stripe phase at higher doping. Within linear spin-wave theory the dynamic 
structure factor is in very good agreement with inelastic neutron scattering data and can account for the 
observed energy dependent anisotropy. 

There is emerging a coherent picture concerning the structure of spin excitations in strongly underdoped 
cuprates (for a recent review cf. e.g. Ref. [1]). In both LSCO (with doping less than 5%)''''' and YBCC 
materials it is now possible to prepare samples with predominantly one twin domain of the orthorhombic 
structure which allows to elucidate the symmetry of the magnetic excitations in the CUO2 planes. Inelastic 
neutron scattering (INS) experiments'"'' ' reveal a strongly anisotropic low energy spin response which breaks 
the four-fold symmetry of a CUO2 plaquette. The latter is of course also broken due to the orthorhombic 
distortion, however, for 'conventional' magnetic excitations this would only induce a small redistribution of 
weight. In cuprates instead, the low energy response is quasi one-dimensional and corresponds to dynamic 
incommensurate spin modulations which are oriented along the diagonals of a CUO2 plaquette in LSCO whereas 
in YBCO the response is oriented along the Cu-0 axis. Upon increasing energy these excitations disperse towards 
the antiferromagnetic wave-vector Q^p which is reached at a characteristic energy Ecross- Above Ecross the magnetic 
fluctuations more or less isotropically disperse outwards and approach the antiferromagnetic (AF) spectrum of 
the undoped compound. On the other hand, a scan in the direction orthogonal to the spin modulation reveals a 
commensurate response up to ficross whereas for larger energies the excitations start to disperse with almost the 
same velocity as in the incommensurate modulation direction. Thus the overall magnetic fluctuation spectrum 
resembles an hourglass parallel to the low energy magnetic spin modulations whereas it is 'Y'-shaped in the 
orthogonal direction. 

Itinerant approaches'" can capture the features of hourglass-shaped excitations, however, it is questionable 
whether these theories are compatible with the strong energy dependent orientational anisotropies described 
above. In fact, the consideration of spin excitation anisotropies in detwinned YBa2Cu306 gs" (i.e. close to optimal 
doping) already requires an anisotropy of —3% in the hopping integrals" along the two orthorhombic directions. 
In contrast, close to the metal- insulator transition orthorhombicity alone cannot account for the strong aniso- 
tropies as also revealed by transport experiments on detwinned samples'*. 

The above experimental results are generally well explained by computations of the dynamical structure factor 
on top of an ordered array of domain walls of the antiferromagnetic order parameter, so-called stripes'^ '". Not 
only the hourglass is well described but also general features as, for example, the recovery of C4 symmetry at high 
energy'\ On the other hand a severe problem with this interpretation for strongly underdoped cuprates is the 
apparent absence of charge order. It has been argued that disorder can 'wash out' the stripe charge correlations 
and corresponding models of disordered (or fluctuating) stripes have in fact been used in order to describe the 
magnetic hourglass spectra of cuprates and cobalt oxides"'^". 

Another kind of theories is based on the formation of spin spirals which also break spin rotational invariance 
similar to stripes, but in the absence of concomitant charge correlations. Early proposals for spirals'" have been 
shown to suffer from an intrinsic instability for homogeneous phases manifested by a negative compressibility^^. 
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On the other hand, holes bound to impurities in the spin glass phase 
may generate a 'dipolar distortion' in the antiferromagnetic spin 
background^"* which in turn can stabilize long-range spiral order^'''^^. 
Such long-range order coexistent with a disordered background is 
reminiscent of liquid crystal phases^*". 

Spin excitations for homogeneous (perfectly ordered) helical mag- 
nets have been computed within the tj-model (see e.g.^^'^") and ana- 
lyzed with regard to the magnetic fluctuations in cuprates^^. 
However, to the best of our knowledge detailed computations of 
the dynamical structure factor for spiral ground states on top of a 
disordered background to compare with experiment have not yet 
been done. 

More recently, we have proposed a state for heavily underdoped 
cuprates which shares characteristics of both the dipolar and the 
stripe scenario. Holes added to the antiferromagnet clump in stripe 
segments ending in fractional vortex- antivort ex (VA) pairs™ resem- 
bling disclination pairs in nematics^"". A single short segment pro- 
duces a dipolar distortion of the AF background similar to Ref [23]. 
Under appropriate conditions a finite density of segments globally 
breaks C4 rotational and inversion symmetry and yields spiral long 
range order, a situation termed ferronematic in Ref [30]. At large 
length scales the spin structure coincides with the state proposed in 
Refs. [24,25], but the microscopic origin of the dipoles, here stem- 
ming from the charge segment nematic order, is completely different. 
In analogy with liquid crystals, at large doping this ferronematic state 
is expected to transform into a smectic phase thus providing a natural 
link with the observed stripe state' around doping 1/8. In fact, the 
ferronematic state has the appealing feature that for parameters 
appropriate for cuprates the resulting magnetic incommensurability 
coincides with that obtained for static stripes so that a smooth cross- 
over from the former to the latter with doping is feasible. 

It has been shown within Gutzwiller variational calculations of 
Hubbard-type models'"'" that VA dipoles correspond to stable 
energy minima for typical cuprate parameters. However, if one is 
interested only in the spin structure of vortex-antivortex pairs it 
turns out'" that the Gutzwiller results can be reproduced within a 
spin-only (Heisenberg) model supplemented by frustrating interac- 
tions in order to stabilize the segments and the induced ferronematic 
order. Such a model may appear too simplistic to compute dynamical 
properties. However, our previous experience with ordered stripes'^ 
shows that it can capture many features of a fully fermionic com- 
putation, which in the present case would be restricted to too small 
sizes to be meaningful. 

Thus in order to study excitations on top of the ferronematic 
ground states our investigations are based on the Heisenberg model 



H-- 



Ny^N^^gV where v denotes the Tilling factor' of a 



■ 0.7 for diag- 



where the spin structure is taken as coplanar and small out-of-plane 
tilts due to the Dzyaloshinsky-Moriya interaction are neglected. The 
segments are modeled as a one-dimensional array of plaquette 
centered magnetic vortices with alternating topological charge (cf 
Fig. la). For a number of N^^g segments the total number of holes is 
given by Nj, 

segment. Variational Gutzwiller calculations™ yield v 
onal segments. 

The following computations are performed for a ratio between 
next-nearest (J') to nearest (/) -neighbor interaction given by /'// = 
— 0.1. These parameters are fixed in the undoped phase by the con- 
dition that the calculated magnon excitations match the measured 
spin-wave dispersion'^. 

Results 

We compute the initial spin structure from a linear superposition of 
randomly positioned vortex-antivortex segments with ferro order as 
in Ref [30] 




(C) 



H M 



(bi ' 


1 






vi 




-J 


© 


— 


1 , 1 


1 





-0.2 -0.1 0 0.1 0.2 
K (r.l.u.) 



Figure 1 | Static properties of the ferronematic spin structure, (a) Detail 
of the staggered ferronematic spin arrangement from a 84 X 84 lattice in 
the xy-plane. Segments are build from plaquette centered vortices 
(squares) and antivortices (circles) which are stabilized by exchange 
couplings ]v between diagonally neighbored spins perpendicular to the 
vortex-antivortex line. Furthermore a nearest-neighbor (/) and next- 
nearest neighbor (/') coupling is implemented in the model, (b) Fit of the 
spin structure factor (LSCO) for doping «/, = 0.05 and segment length 
= 4 perpendicular to the direction of segments. Computations have been 
done on 84 X 84 lattices and we average over 10 segment configurations. 
The experimental resolution (cf. Fig. 4 of Ref [2]) has been taken into 
account by convoluting with a gaussian. Data by courtesy of S. Wakimoto. 
(c) Sketch of the Brillouin zone geometry used in the present paper. 
The JC-direction corresponds to the direction of the disordered spiral 
modulation and the shaded square indicates the momentum range for the 
intensity plot of Fig. 3. 



S'ir) = So exp(!Qr)cos a>(r) 
y(r)=Soexp(!Qr)sin(l)(r) 
S^(r) = 0. 

Here Q = (tt, tt) is the AF wave-vector and 
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arctan 
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is the angular spin distortion due to the Nj^g segments. In Eq. (3) the 
sum over i{n) is due to sites belonging to segment n and ( — 1)"'"* = 
1(— 1) corresponds to a plaquette centered magnetic (anti) vortex. 
Then the classical energy of Eq. (1) is minimized where the VA 
structure is stabilized by an additional coupling ]v between diag- 
onally neighbored spins across the segments (cf Fig. la). A fit of 
the VA structure to the result of a Gutzwiller variational computation 
yields Jyl] = 8 (cf Methods section). In order to account for fluctua- 
tions beyond the variational result we take a smaller value Jy/J = 1.5 
which induces the same winding number for the ferronematic spiral 
as for Jy/J = 8 (cf supplementary information) and reproduces the 
low doping magnetic incommensurability of LSCO. Thus this para- 
meter is fixed at the static level and there are no free parameters in the 
following dynamical computations. 

For doping of n/, = 0.05 Fig. lb demonstrates that this approach 
can reproduce the static spin structure factor of LSCO measured with 
elastic neutron scattering. Note that the present fit with plaquette 
centered vortices yields shorter segments with N^,^ = 4 as compared 
to our previous investigations'" with site-centered vortex-antivortex 
pairs. 
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Figure 2 | Neutron scattering intensity for the ferronematic at doping nj, = 0.05. The magnetic neutron-scattering cross section has been integrated over 
energy ranges (a) 0-0.06 /, (b) 0.13-0.2 /, (c) 0.27-0.34 /, (d) 0.68-0.75 /. The momentum range corresponds to the shaded area indicated in Fig. Ic. The 
value of the perpendicular momentum q^has been calculated from energy-momentum conservation for an incident neutron beam energy with £; = 0.3 / 
(panels a,b) and £, = 0.8 / (panels c,d). Average over 6 configurations. 



In order to compute excitations on top of the ferronematic ground 
state we rotate the spin quantization axis to the direction of the local 
mean-field magnetization and apply linear spin-wave theory as 
described in Methods. 

It is then straightforward to obtain the spin susceptibilities 

; Sf)) from the equations of motion and to compute 



the Fourier transformed 



config. 



(4) 



In practice we evaluate y^\oj) for a randomly generated configura- 
tion of segments and then take the average over typically 10-20 
configurations as indicated by (. . ■) config- 

The spins of the ferronematic state can fluctuate both in- and out- 
of the Cu02 plane. Both types of excitations contribute to the mag- 
netic neutron- scattering cross section measured by INS which is 
proportional to 

dQ.dE ' 



Here x°'''(q> only depends on the in-plane components of q while 
the perpendicular component enters the cross section only 
through the prefactors inside the sum. In the extreme case that the 
is much larger than the in-plane components (i.e. for the case that 
the neutron beam is almost perpendicular to the Cu02 plane and 
high energies) the response is dominated by the in-plane response 
y}^{oS)^^. In practice, q^ depends on energy and one always has a 
contribution from perpendicular, in-plane and off-diagonal suscept- 
ibilities which we have all included. In our computations we assumed 
an incoming neutron energy £, and a scattering geometry as in typ- 
ical neutron scattering experiments (i.e. q^ and qy in the first Brillouin 
zone). 

Fig. 2 shows momentum cuts of the neutron scattering intensity 
within different energy ranges which have been obtained from the 
same ferronematic configurations for which the static response in 
Fig. lb has been evaluated. Momentum cuts along the H- and 



K-direction for the same energies are reported in Fig. 3. At low 
energies (panel (a) in Figs. 2, 3) the response is dominated by two 
incommensurate peaks along the orthorhombic K-direction due to 
the spiral modulation perpendicular to the segments and coinciding 
in momentum space with the static peaks in Fig. lb. The period of the 
modulation is obtained as a = 27c/|gx| = lAUcrtho corresponding to 
an incommensurability of fi"""" = 1/14 = 0.07143. Note that for 
diagonal stripes with v holes per site the incommensurability would 
be given by e°"'"' = n/,/v thus yielding the same value for ni, = 0.05 
and V = 0.7. 

Upon increasing energy the reponse first keeps the elongated 
shape along the K-direction (panel b) but becomes more isotropic 
at higher energies (panel c). From the momentum cuts in Fig. 3 the 
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Figure 3 | Momentum cuts parallel (//-direction, red squares) and 
perpendicular to the segments (if-direction, black circles) for the 
intensity distributions shown in Fig. 2. 
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Figure 4 | Intensity plot of the cross section (oSq(co) [Eq. (5)] for a 
momentum cut perpendicular (panel a) and parallel (panel b) to the 
segments. The value of the perpendicular momentum has been 
calculated from energy-momentum conservation for an incident neutron 
beam energy vnth £, = /. The solid lines report the magnon dispersion of 
the undoped antiferromagnet. Average over 6 configurations. 

dispersion of excitations along the K-direction toward the antiferro- 
magnetic wave- vector can be clearly deduced. Panels b,c of Fig. 3 also 
reveal a small shoulder in the momentum structure along K which is 
due to an outwards dispersing branch. Similar to the case of stripes"" 
and enhanced by disorder the corresponding intensity is strongly 
reduced with respect to the inwards dispersing excitations. 

At stiU higher energies the magnetic scattering intensity acquires a 
quasi two-dimensional ring-shaped structure (Fig. 2, panel d). The 
overall features of the intensity distribution are rather similar to data 
from INS experiments on LSCO (cf Fig. I in Ref. [9]), however, our 
calculations reveal an anisotropic intensity distribution with variable 
weight ratio between the H- and K-directions at different energy 
ranges, (cf. panel (c) of Fig. 3). This feature resembles an analogous 
rotation of weight observed in the magnetic scattering from static 
stripes in Lai 875Bao.i25Cu04'''', whereas the data from strongly 
underdoped LSCO seem to indicate a more isotropic distribution 
of weight. We should remark that for the case of stripes we observe 
that a computation in a spin only model (as here) yields a more 
anisotropic response than a more accurate computation in a fermion 
model'"". Thus we attribute this small difference in anisotropy to the 
simplicity of the model. Unfortunately a fermion model could not be 
used due to the large size required to take into account disorder with 
high momentum resolution. 

The magnetic dispersions along K and H directions can be 
deduced from Fig. 4 which shows the corresponding cuts of the 
magnetic structure factor Eq. (5) weighted by energy. Panel (a) of 
Fig. 4 clearly reveals the basic features of the hourglass spectrum 
although the resolution is limited due to the large periodicity of the 
ferronematic spiral as compared to the size of the lattice which can be 
dealed with. The large intensity feature corresponds to the inwards 
dispersing magnetic excitations which merge at Q^p and subse- 
quently disperse outwards again with reduced intensity. At high 
energies the magnetic fluctuations approach the magnon dispersion 
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Figure 5 | Energy dependence of the magnetic cross section Sq=Q„(co) 
(cf. Eq. (5)) at the antiferromagnetic wave-vector. The main panel reports 
Sq=Q^,(co) for doping x = 0.04 (dashed) together with experimental data 
extracted from Fig. 1 of Ref [7] . Inset: Comparison of the crossing energy 
with experimental data for LSCO (solid circles) and LBCO (open circle) 
from Ref. [8] . For the conversion into energy units we use / = 112 meV^ 
and a frequency independent renormalization factor = 1.18''^. 

of the undoped antiferromagnet which is indicated by the solid line 
in Fig. 4. The Y-shaped structure of the spin excitations in the parallel 
direction to the segments is reported in panel (b) of Fig. 4. Here the 
low energy magnetic excitations are confined to Qj^p but also start to 
follow those of the undoped AF above the crossing energy E„oss- 

The location of Ecross is difficult to determine from Fig. 4 due to the 
limited resolution. Instead we deduce this scale from the energy 
dependence of the magnetic structure factor for the antiferromag- 
netic wave vector Qj^p which is shown in Fig. 5 for doping x = 0.04. 
Since the ferronematic spin spiral is strongly disordered the static 
structure factor also displays small but finite weight at Qaf- This 
induces a 'quasi-elastic' peak at oj ^ 0 in the frequency dependent 
magnetic cross section, Sq=Qj^((a) in addition to the higher energy 
feature at E^^^^ ~ 0.14 /. We can convert the 'J-scale' into energy 
units by taking / = 112 meV'^ and a frequency independent 
renormalization factor Z^= 1.18'^ which yields good agreement with 
the magnon dispersion of the undoped AF. Based on this conver- 
sion we also show in Fig. 5 experimental Sq = Q^f(aj) data for 
Lai 96SrOo.o4Cu04 which we have extracted from Fig. I of Ref [7]. 
Interestingly also these data support an increase of Sq=Qjp((») for co 
^ 0 in agreement with the persistence of AF domains in an otherwise 
incommensurate ordered spin system. Finally, the inset to Fig. 5 
reports Ecross for X = 0.04 and x = 0.05 which fits rather well into 
the experimentally determined" evolution of £cross with doping. In 
contrast, stripe calculations in the strongly underdoped regime usu- 
ally significantly overestimate Ecross" although this may also be cured 
by disorder". 

Discussion 

We have shown that the spectrum of spin excitations computed on 
top of a ferronematic ground state is in qualitative and quantitative 
agreement with inelastic neutron scattering experiments on cuprate 
superconductors in the spin glass phase. Since the ferronematic also 
can be thought of arising from melted stripes our theory provides an 
explanation for the continuous evolution of the hourglass-shaped 
magnetic excitations from low to optimal doping. It should be 
mentioned that preliminary Monte-Carlo simulations also seem to 
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Figure 6 | Comparison of spin textures obtained from the Hubbard model within the GA (a) and the Heisenberg model (b) The GA solution shown in 
panel (a) has been obtained for 4 holes in a 14 x 14 system. Parameters: U/t = 8, t'/t = —0.2. The spin structure of panel (b) results from the 
minimization of the Heisenberg model. Here the 4 diagonal exchange couplings which are indicated by bars are set to J^/J = 8 and aU other nearest- 
neighbor couplings are /'// = —0.1. Panels (c), (d) compare the angular rotation of the spins ( = atan(Syl S^)) along the scans indicated in panel (b). 



support smectic charge correlations which are stabilized by long- 
range Coulomb interactions. While such correlations between dis- 
ordered stripes can also induce an incommensurate magnetic 
response''^, the magnetic fluctuations of the ferronematic spiral are 
due to the disordered magnetic background which is induced by the 
vortex and antivortex dipoles of the charge segments. 

Experimentally, the ferronematic state is supported by the obser- 
vation of a low temperature ferroelectric phase in the spin glass phase 
of cuprates"". In fact, the ferronematic state breaks inversion sym- 
metry so that on general grounds'' a ferroelectric distortion becomes 
possible e.g. via the inverse Dzyaloshinskii-Moriya mechanism""*. 
The 'smoking gun' experiment would then be to test whether the 
onset temperature of incommensurate spin fluctuations coincides 
with that of the ferroelectric as predicted by our model. 



(indicated by circles) are arranged in 2 X 2 plaquettes aligned along the diagonal 
direction. The center of the plaquettes alternately contains a magnetic vortex and 
antivortex, respectively, inducing a phase change of the AF order across the segment. 
This feature can be modeled within the Heisenberg model by diagonal couplings 
indicated by solid bars in panel (b) of Fig. 6. Panels (c,d) compare the angular spin 
twist — atan{SyfSx) between GA computation and minimization of the Heisenberg 
model energy along the two scans indicated in panel (b). We find that for Jy/] — 8 
there is excellent agreement of the spin far field {which is determined by the dipole 
strength of the segment) between both approaches. On the other hand, one may argue 
that the variational GA overestimates the localization of hole charges and that 
including corrections beyond the G A would lead to actually smaller values of /,,//. The 
phase change across the segment is then reduced for smaller /v//as can be seen from 
panels {c,d). However, as long as Jy/J £ 1.5 we find that the static and dynamic spin 
structure factor only weakly depends on ]J] as demonstrated in the supplementary 
material where we show results for the GA fitted value 7^/7—8. In the 'Results' section 
instead we have taken a smaller value 7,-/7 — 1.5 in order to include fluctuations 
beyond the GA variational approach. 



Methods 

Parameters for the effective Heisenberg model. In order to estimate the effective 
coupling /,, across tlie plaquette centered magnetic vortices and antivortices (cf. Fig. 1 ) 
we perform a variational calculation based on the single-band Hubbard model 

including nearest- neighbor (t) and next-nearest neighbor {t') hopping. Ground state 
properties are obtained from the Gutzwiller variational wave-function ( Tc; \—Pg\^o) 
where Pc- — H, [l — (1 ~^i)"i,T"(u] partially projects out doubly occupied sites from 
the Slater determinant |Oo). We evaluate ^ <*i'G|^|*I'G>/(*i'G|^G> in the limit 

of infinite dimensions'^ corresponding to the so-called Gutzwiller approximation 
(GA) and minimize with respect to the variational parameters gf. We set U/t ~ 8 and 
t'/t = —0.2 which within a time-dependent extension of the GA reproduces the spin- 
excitations of undoped LSCO^^'^^. 

Panel (a) of Fig. 6 reports the spin/charge texture obtained in this way for 4 holes 
(192 particles) on a 14 X 14 lattice with periodic boundary conditions. Hole charges 



Evaluation of spin excitations. Spin excitations on top of the ferronematic ground 
state are obtained by applying a rotation to the spins 

Sj^ - cos(0;) -\- ^ sin(0,-) (7) 
S^'-S;'cos(0i)-FSf sin(0,) (8) 



Si-~S- (9) 

so that the spin quantization axis points in the direction of the local mean-field 
magnetization^". 

For this effective ferromagnetic state one can apply the Holstein-Primakoff 
transformation: 
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2i\ ' 

The linear-spin wave theory Hamiltonian reads. 

1 



with 



•i 'i 



,^j=-I,j[cos{@,-&^)- 



By=-hj[cos{@i-&j] 



-1] 



(10) 



(11) 



(12) 



(13) 



where 0, is the angle the spin i forms with the x-axis. Here we used the fact that the 
configurations are coplanar. The absence of unstable spin-wave modes in our analysis 
supports this assumption. 
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